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We explore the prospects for indirect detection of neutralino dark matter in supersymmetric mod- 
els with an extended Higgs sector (NMSSM). We compute, for the first time, one- loop amplitudes 
for NMSSM neutralino pair annihilation into two photons and two gluons, and point out that extra 
diagrams (with respect to the MSSM), featuring a potentially light CP-odd Higgs boson exchange, 
can strongly enhance these radiative modes. Expected signals in neutrino telescopes due to the 
annihilation of relic neutralinos in the Sun and in the Earth are evaluated, as well as the prospects 
of detection of a neutralino annihilation signal in space-based gamma-ray, antiproton and positron 
search experiments, and at low-energy antideuteron searches. We find that in the low mass regime 
the signals from capture in the Earth are enhanced compared to the MSSM, and that NMSSM 
neutralinos have a remote possibility of affecting solar dynamics. Also, antimatter experiments are 
an excellent probe of galactic NMSSM dark matter. We also find enhanced two photon decay modes 
that make the possibility of the detection of a monochromatic gamma-ray line within the NMSSM 
more promising than in the MSSM. 



I. INTRODUCTION 



Numerous theoretical and phenomenological motivations exist for a Minimal Supersymmetric extension of the 
Standard Model (MSSM). At the same time one of the attractive by-products of low energy supersymmetry is the 
natural occurrence in the particle content of the theory of a stable weakly interacting massive particle, the lightest 
neutralino, which could be the microscopic constituent of the as yet unobserved galactic halo dark matter. Another 
strong motivation comes from the SM hierarchy problem, originating from the large fine-tuning required by the 
stability of the electroweak scale to radiative corrections, originating from the large number of orders of magnitude 
occurring between the GUT, or the Planck, scale, and the electroweak scale itself. 

Although very appealing, the MSSM has been challenged by various pieces of experimental information, and by 
some arguments of more theoretical nature. Among these, the LEP-II limit on the mass of the lightest CP-even 
Higgs jl|, the constraints on the masses of supersymmetric (Susy) charged or colored particles from direct searches 
at LEP and at the Tevatron 2j, and the so-called problem, i.e. the fundamental reason why the SuSY Higgsino 
mass term fi appearing in the MSSM superpotential lies at some scale near the electroweak scale rather than at some 
much higher scale. 

The addition of a new gauge singlet chiral multiplct, S, to the particle content of the MSSM can provide an elegant 
solution to the mentioned /i problem of the MSSM 3]. The so-called Next to Minimal Supersymmetric Standard 
Model (NMSSM) is an example of one such minimal extension that also alleviates the little fine tuning problem of 
the MSSM, arising from the non-detection of a neutral CP-even Higgs at LEP-II 

A further motivation to go beyond the MSSM comes from Electro- Weak Baryogenesis (EWB), i.e. the possibil- 
ity that the baryon asymmetry of the Universe originated through electro-weak physics at the electro-weak phase- 
transition in the Early Universe. Although still a viable scenario within the MSSM EWB generically requires the 
Higgs mass to be in the narrow mass range above the current LEP-II limits and below ~ 120 GeV, a rather unnatural 
mass splitting between the right-handed and the left-handed stops (the first one required to lie below the top quark 
mass, and the other in the multi-TeV range), CP violation at levels sometimes at odds with electric dipole moment 
experimental results, and, generically, a very heavy sfermion sector |^. In contrast, the NMSSM provides extra 
triscalar Higgs couplings which hugely facilitate the occurrence of a more strongly first-order EW phase transition, 
and extra CP violating sources, relaxing most of the above mentioned requirements in the context of the MSSM (Te,^. 

One of the chief remaining cosmological issues associated with the NMSSM, the cosmological domain wall prob- 
lem '9^, caused by the discrete Z3 symmetry of the NMSSM, can be circumvented by introducing non-renormalizable 
Planck-suppressed operators [Tol |. 

The Higgs sector of the NMSSM contains three CP-even and two CP-odd scalars, which are mixtures of MSSM-like 
Higgses and singlets. Also, the neutralino sector contains five mass-eigenstates, instead of the four in the MSSM, 
each of which has, in addition to the four MSSM components, a singlino component, the latter being the fermionic 
partner of the extra singlet scalars. The extended Higgs and neutralino sectors weaken the mass bounds for both the 
Higgs bosons and the neutralinos. Very light neutralinos and Higgs bosons, even in the few GeV range, are in fact not 
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excluded in the NMSSM (see also the particle spectrum with the dominant 1-loop and 2- loop corrections 

to the Higgs sector is available via the numerical code Nmhdecay jl^). 

The cosmology of Dark Matter singlinos has been addressed long ago in Ref. ) while the LSP relic abundance in 
the NMSSM was first calculated in jl^. Constraints from electroweak symmetry breaking and GUT scale universality 
were added in (1^. More recently, the computer code micrOMEGAs ilTj], has been extended to allow for the relic 
density calculation in the NMSSM pl|^ 

The implications for the direct detection of NMSSM neutralinos were first studied in l20l|, where light neutralinos 
3 GeV) with acceptable relic abundance and sufficiently large expected event rates for direct detection with a 
^^Ge-detector were found in different domains of the parameter space, when the gaugino unification relation was 
relaxed. In |2l|. the theoretical predictions for the spin-independent neutralino-proton cross section, (Tj^o_p, were 
reevaluated, and all available experimental constraints from LEP on the parameter space were taken into account. 
Values within reach of present dark matter detectors were obtained in regions with very light Higgses, m^o ^ 70GeV, 
with a significant singlet contribution. The lightest neutralino, in those regions, features a large singlino-Higgsino 
composition, and a a mass in the range 50GeV < m^o < lOOGeV. More recently, NMSSM neutralinos as light as 
lOOMeV ^ m^o < 20GeV, satisfying accelerator constraints and with the right relic density, have been shown to occur 
in l^^, where it was argued that the NMSSM can, moreover, provide neutralinos in the mass range that would be 
required to reconcile the DAM A claim of discovery with the limits placed by CDMS. 

So far, theoretical studies have not addressed the possibility that NMSSM neutralinos making up the galactic dark 
matter can manifest themselves indirectly. For instance, pair annihilations of neutralinos in the galactic halo can 
produce sizable amounts of antimatter, which current and forthcoming space-based antimatter search experiments 
can possibly detect; neutralinos trapped in the Sun or in the Earth 23, 24] can give rise to a coherent flux of energetic 
neutrinos from the center of the Sun or of the Earth; pair annihilation of neutralinos, either in nearby large-dark- 
matter-density sites, or from the cumulative effect of annihilations outside the Galaxy, can produce gamma-ray fluxes 
at a level detectable by GLAST or by ground-based air Cherenkov Telescopes; and last, but not least, the exciting 
possibility of peculiar gamma-ray spectral features, like a sharp monochromatic peak at Ej ~ (where x indicates 
the lightest neutralino, assumed to be the lightest supersymmetric particle), from loop- induced xx ~* 77 processes, 
can also be, in principle, very promising. 

A first motivation for looking into indirect dark matter detection within the NMSSM comes from the possibility that 
thermally produced neutralinos, in this context, can be very light. Since the pair annihilation rate of thermal relics 
is roughly fixed by requiring that the thermal neutralino abundance coincides with the cold dark matter abundance 
inferred by astrophysical observations [2^, the indirect detection rates generically scale as 1/m^: light neutralinos 
are therefore expected to give significantly enhanced rates with respect to the standard case. 

A more technical point has provided us with a second motivation to look into indirect detection prospects for 
NMSSM neutralinos: loop induced pair annihilations of neutralinos into two photons or two gluons (respectively 
contributing to the mentioned monochromatic gamma-ray line and to, e.g., antimatter fluxes) are predicted to be 
increased, within the NMSSM, by the presence of extra diagrams mediated by the (potentially light) extra CP- 
odd Higgs boson. We therefore extend here, for the first time, the MSSM results for these loop-induced neutralino 
pair-annihilation amplitudes [2^ to the NMSSM. 

Our results suggest several signatures, including muons resulting from neutralino annihilation in the Earth, and 
antiparticle and gamma ray production from neutralino annihilation in the galaxy, where the NMSSM produces signals 
that are enhanced compared to those predicted in the MSSM. 

The outline of this article is as follows: we first introduce the theoretical framework and set our notation in Sec.lTTI 
we devote Sec. Ill Al to a discussion of the viable NMSSM parameter space. Sec. IIIII contains our central results 
on indirect NMSSM neutralino dark matter detection, while the appendices provide the reader with details on the 
relevant NMSSM neutralino pair- annihilation amplitudes and on neutralino- nuclcon scattering cross sections. 

II. THE NMSSM 

We hereby describe the Lagrangian of the NMSSM. Our notation follows that of the code Nmhdecay 0|, which 
we have used to explore the NMSSM parameter space^. 



^ The pheno meno logy of the hghtest neutralino in a different extension of the MSSM, the Left-Right SuSY model, has been recently 
surveyed in Il9l . 

^ Note that the Higgs states Hu, are usually denoted in the MSSM by H2 and Hi. As shown in the Appendix, some indices in both 
the neutralino and the Higgs mass matrices need to be switched accordingly to make contact with the corresponding MSSM expressions. 
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Apart from the usual quark and lepton Yukawa couplings, the scale invariant superpotential is^ 

A SHuHd + 2 



(1) 



depending on two dimensionless couplings A, k beyond the MSSM, and the associated trilinear soft-SuSY-breaking 
terms 



(2) 



The two other input parameters, tan/3 = {Hu) / {Hd) and ^cS = A(S'), along with Mz, determine the three SuSY 
breaking masses squared for Hu, Hd and S through the three minimization equations of the scalar potential. Note 
that an effective yit-term is generated from the first term in Eq. ^ for a non-zero value of the VEV (5). With the 
sign conventions of 'l3] for the fields, A and tan/3 are positive, while k, A\, A^ and /^cff can have either sign. 

Assuming CP conservation in the Higgs sector, there is no mixing between CP-even and CP-odd Higgses. More 
concretely, for VEVs /i„ = hd = {Hd) and s = (S) such that 
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(3) 



the CP -even mass matrix in the basis S" "'^ — {Hur, HdR, Sb) is rendered diagonal by an orthogonal matrix 5"^ . 
One, thus, obtains 3 CP-even mass eigenstates hi — SijSj""^'^, with increasing masses mh.. The bare CP-odd states 
pbare _ (^ff^j^ff^j^ gj^ are related to the physical CP-odd states a^, i = 1, 2, and the massless Goldstone mode 03 = G 
by fli = Py Pj""''^, where oi and 02 are ordered with increasing mass. Details of the bare mass matrices in terms of 
the NMSSM parameters can be found in Ts'l. 

With fixed parameters of the Higgs sector, the masses and mixing of the neutralinos are determined by two additional 
parameters: the masses Mi and M2 of the U{1)y gaugino, Ai, and the neutral SU{2) gaugino, A2. In the basis 
■0'^ = {—iXi,—iX2,ip^,ipdTips) the symmetric mass matrix A4q of the neutralinos 



(4) 



has the form 
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This matrix can be diagonalized by a real orthogonal matrix, A^ij, obtaining 5 eigenstates, 
but not necessary positive masses, m^o, ordered in increasing absolute value of the mass^. 



Nijip'j, with real, 



A. Light neutralino dark matter: parameter space 

Even though our study of the indirect detection of NMSSM-like neutralinos is completely general, we choose to 
focus on light neutralinos, m^o ^ lOOGeV, since the differences with the case of the MSSM will be more acute in this 
case. 

We have performed a scan of the parameter space with the program Nmhdecay. For each point, after computing 
the masses and couplings of all physical states in the Higgs, chargino and neutralino sectors, Nmhdecay checks for 
the absence of Landau singularities below the GUT scale for A, k and the Yukawa couplings ht and h^. This translates 
into A < .75, k < .65, and 1.7 < tan/3 < 54 Nmhdecay also checks for the absence of an unphysical global 



^ Hatted capital letters denote superfields, and unhatted ones their scalar components. 

* The matrix jFJ can also be diagonalized using a complex Nij. In that case the mass eigenstates would b e re al and positive. These two 
choices result in diflferent signs of certain Feynman rules, as pointed out in the Appendix (for details see l27l V 
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minimum of the scalar potential with vanishing Higgs VEVs^. The program also makes sure that Higgs and squark 
masses are positive, thus avoiding, in particular, charge breaking minima. 

Finally, the available experimental constraints from LEP are imposed, including unconventional channels relevant 
for the NMSSM Higgs sector, bounds on the invisible Z width (for light ncutralinos) and limits on chargino and 
neutralino pair production. 

As remarked in ;21j, there are other experimental bounds that might put constraints on the parameter space. 
Rare B-meson decays, sensitive to physics beyond the Standard Model like supersymmetry, have been studied for the 
NMSSM in the large tan/3 regime 28J. However, the transitions 6 ^ 57 and Bg are both flavor changing 

while the contributions from a light can be suppressed by making the appropriate squark or slepton heavy |22| . 

Additional constraints, that apply when the x° and oi are light, were considered in [22|. The conclusion is that 
bounds on the magnetic moment of the muon can only be violated in extreme models, while rare kaon decays rule 
out some models with extremely light a\. 

On the other hand, decays of the vector resonances J and T might be important for models with a light 
and/or a\ (we follow here the discussion in H^]). In some of our models, the decay V aij, where V stands for J/'i' 
or T, is indeed possible^. The width, relative to the muon decay channel is at leading order p9l |: 



r(F->7ai) ^ Gprn l n? 



\-^\X\ (6) 



where X — tan fiP'^i for the T decay, X — cotan/3P{]^ for the J/^ and the gives the piece of a\ that would be the 
MSSM pseudoscalar if the singlet were not present^. 

The ratio in Eq. © is generally less that 4 x 10^^ (0.006) for J/^ (T) decays, below the CLEO measurement 
of T ^ invisible + 7 [s^l, although the highest branching ratios could be discovered with new upcoming data or 
reanalyzing the existing CLEO data On the other hand, lepton universality tests in T decays by high- luminosity 
B factories could detect a CP-odd Higgs in the mass range 5 GeV ^ rua^ ^ mx, which would be otherwise very hard 
to discover by just looking at the 017 channel [sj. 

To generate our models, we scan, at random, the NMSSM parameter space in the region: 

< A < 0.75 

-0.65 < At < 0.65 

1.7 < tan/5 < 54 

80GeV < [I < 500GeV 

-500 < Ax,A^ < 500. (7) 

The gaugino masses were also randomly chosen within the bounds OGeV < AIi < 100 GeV, Mi < M2 < 500 GeV and 
300 GeV < M3 < 1000 GeV. The soft sfermion masses were set to Afa, and the sfermion trilinear terms were varied 
within ±1.5M3. 

The models that passed the phenomenological constraints imposed by Nmhdecay were fed into micrOMEGAs 
to calculate the Xi relic density, taking into account all possible annihilation and coannihilation channels. We kept 
as viable those models that fell within the 2 — a region for the Cold Dark Matter abundance inferred by the WMAP 
team for a ACDM cosmology psj . 

We show in Fig. ^ the region of the NMSSM parameter space that satisfies the constraints discussed above. We 
have classified an NMSSM-like neutralino as bino-like if Nfi > 0.9, singlino-like if fulfilling the condition > 0.9, 
otherwise we indicate the neutralino "mixed" . 

As can be seen from Fig. ^ the neutralino is mostly a singlino when k and A are small. We can understand this 
feature by realizing that the upper 4x4 block in the neutralino mass matrix, Eq. ||SJ), corresponds to the MSSM. 
From the lower 3x3 block it can be appreciated that the singlino decouples from the MSSM part when [s^: 

2\ks\, Xv < Ml, M2, \fi\. (8) 

It must be stressed that our Fig. ^ shows only models which give an acceptable relic density. This might explain 
the absence of singlino-like ncutralinos at moderate A ~ 0.3. 



^ Tree level restrictions in parameter space leading to valid minima are discussed in |21| . 

® Our scan of the NMSSM parameter space does not yield Xi light enough to make the d ecay V —> 7X1X1 kinematically allowed. 

The definition of in terms of the CP-odd Higgs mixing matrix, Pij , can be found in Il3l . 
* We find one model in our scan, with = 3.01 GeV and tan/3 = 48.24, yielding a ratio for the T decay large enough to be already 

excluded by CLEO. 
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FIG. 1: The allowed region of (A — k) parameter space (upper-left); the m-o vs. plane (upper-right)- the decay ai — > XiXi 
is allowed for those models below the line; the (tan /3 — A) viable parameter space (lower-left); as a function of the lightest 
CP-odd Higgs mass (lower-right) 



A CP-odd Higgs so light that the decay ai XiXi is possible, is generally not viable for light singlino-hke neu- 
tralinos. This makes them cosmologically disfavored, since this resonant decay is required to enhance the annihilation 
cross-section and obtain the correct relic density. A nearly complete mass degeneracy between Xj* and either the 
next-to- lightest neutrahno, or the lightest chargino, suppressing the LSP final relic abundance through large 
coannihilation effects, usually occurs for the viable light singlino models 17]. 

We can see from Fig. ^ that singlino models at large A feature small tan/3 < 5, since large values of tan/3 induce 
sizable singlino mixing. We also expect models with moderate tan (3, for which annihilation through a Higgs resonance 
is marginal in the MSSM, to be pecuhar of the NMSSM setup. 

Finally, the lower-right panel in Fig. ^ shows how a light CP-odd Higgs boson, oi, appears when — + 0. We 
stress that this regime, when the U{l)ii symmetry approximately holds, is well motivated in the context of gaugino 
mediated SuSY breaking Q where A^, is only generated at the two loop level. 



III. NMSSM DARK MATTER INDIRECT DETECTION 



Neutralinos, being weakly interacting and electrically neutral particles are very difficult to observe in collider 
experiments directly. If they make up a sizable fraction of the galactic halo dark matter, however, other methods of 
detection become feasible [33.l35| . 

Monitoring the energy deposited as they scatter off nuclei in detectors falls into the realm of direct detection methods. 
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A group of experiments is actively exploring this path, although, as already mentioned, their sensitivity decreases 
for neutralinos below 100 GeV^. The prospects for direct detection of NMSSM neutralinos have already been 

discussed in the literature [23, lU |23| • 

We focus here on the possibility that dark matter neutralinos can be detected by looking at products of their pair 
annihilation. Chief among them are neutrinos, photons and antiparticles |34. i35t i36t. i37i | . 

Neutrino fluxes from neutralino annihilations are searched for in underground neutrino telescopes. Present facilities 
such as Super-Kamiokande and MACRO, with low energy thresholds, ^ 1 GeV, are particularly useful to constrain 
the light NMSSM dark matter particles that we consider. Some of the planned facilities (e.g. AMANDA, ICECUBE) 
are geared to detecting neutrinos above 100 GeV, and we don't consider them here. ANTARES, on the other 
hand, promises to improve the sensitivity to moderately energetic neutrinos by an order of magnitude and will be of 
importance for our discussion. 

Gamma-rays are also produced in neutralino annihilations. They can be detected by earth-based Cherenkov tele- 
scopes (MAGIC, HESS, VERITAS, . . . ) or in space-borne facilities (EGRET, GLAST, AMS), although, only the 
latter have the ability to observe the low energy photons from light neutralinos such as those considered here. 

Other satellites, like PAMELA, GAPS and AMS, will measure the flux of antiparticles and antimatter nuclei. 

We study below the signatures of light NMSSM dark matter particles in neutrino telescopes, gamma-ray satellites 
and antimatter detectors. We also touch upon the effects on the Sun caused by neutralino energy transport. 

The main ingredient for indirect detection prediction are the different annihilation modes of neutralinos. Since 
dark matter in the halo moves at non-relativistic velocities, v ~ lO^'^c, only the channels with a CP-odd final state 
can occur. The branching ratios for the relevant tree level processes are reviewed in App. ^ together with the most 
important one loop channels. 

Neutrino fluxes from neutralino annihilations are enhanced in the direction of the center of the Sun or of the Earth. 
The abundance of neutralinos trapped within these objects depends on the scattering cross sections of neutralinos 
with nuclei that can be found in App. IbI 

We start our discussion by studying the information that can be gained from the observation of neutrino fluxes. 



A. Neutrino Fluxes from Neutralino annihilations in the Earth and in the Sun 

The observation of energetic neutrinos from annihilation of neutralinos in the Sun |23[|38j| and/or the Earth [2j 
is a promising method for indirect detection of neutralino dark matter (see e.g. [s^. l35| for a review). 

Neutralinos making up the dark matter in the halo of the Galaxy have a small but finite probability of elastically 
scattering from a nucleus in a given body (the Sun or the Earth). In doing so, neutralinos might be left with a velocity 
smaller than the escape velocity and, thus, become gravitationally bound to the body. The captured neutralinos settle 
to the core of the body, via additional scatterings from nuclei in the body, and eventually annihilate with one another. 

The pair annihilation of the accumulated neutralinos generates, via decay of the particles produced in the various 
annihilation final states, high energy neutrinos with a differential flux given by: 



BRf-^ (9) 



f 

Here d is the distance of the detector from the Sun or the center of the Earth, Tann is the annihilation rate of the 
neutralinos, BRf is their branching ratio into the final state xx ~^ f ^^'^ dNf/dE is the neutrino spectrum from the 
decay of the particles in the final state /. 

Since neutralinos inside the Sun or the Earth are highly non-relativistic, their annihilations occur almost at rest. 
The branching ratios of the different annihilation channels are discussed in App. El 

A light neutralino can only annihilate into the light quarks and lepton pairs, which, after decay, give rise to a fairly 
soft neutralino spectrum. A more massive x can lead to W'^W^ , ZZ and heavier quark pairs, which typically produce 
a harder differential neutralino flux. Apart from these fundamental channels, neutralino annihilations can produce 
Higgs bosons or mixed Higgs/gauge boson final states. The Higgs bosons will, in turn, decay to other Higgses or to 
one of the "fundamental" channels |40| . In our calculations, wc have taken into account the fact that the number of 
final states containing Higgs bosons is increased in the NMSSM due to the extra CP-even and CP-odd states, and 
02, compared to the MSSM. 



In this section we denote the lightest neutralino, Xii simply by x- 
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1. Annihilation rate in the Sun and in the Earth 



Neutralinos accumulate in the Sun or the Earth by capture from the halo of the Galaxy, and are depleted by 
annihilation and by evaporation. The evolution equation for the number of neutralinos, N, in the Sun or the Earth 
is given by: 

^ = C - CaN'^ - CeN (10) 

where C is the rate of accretion onto the body, the second term is twice the annihilation rate and the last term 
accounts for neutralino evaporation. 

Evaporation has been shown to be important only for neutralinos lighter than 3 — 5 GeV [4ll li^ . The lightest 
neutralino that we consider in this paper is on the upper range, to^ ^ 5 GeV, so we can safely neglect the last term 
in Eq. ((TU|l. 

We then solve Eq. (|10|l for N , and obtain the annihilation rate at any given time: 

^ann = ^tauh^ (</t^) (11) 



where ta — 1 / V CCa is the time scale for capture and annihilation equilibrium to occur. We will be interested in the 
value of Tann today, for t = ~ 1.5 x 10^'' s. 

The annihilation rate per effective volume, Ca is given by: 

Ca = {cJAv) (12) 

and Vj — [3TOp;T/(2jm^p)] are the effective volumes for the Sun [Vj ^ 6.6 x 10^*(j'm;^_io)~^/^ cm'^) or the Earth 
{Vj ~ 2.3 X W^{jm^,w)~^'^ cm3). 

The total annihilation rate, {gav), is calculated with all the contributions at tree level, with the inclusion of the of 
the two gluon channel discussed in App. ^ 

The accretion rate in the Sun was first calculated in [2^ and , and for the Earth in [23. Is^ . More detailed 
evaluations can be found in [4^ • The results depend on the velocity dispersion in the halo, the velocity of the Sun 
with respect to the halo, the local density of dark matter and the composition of the Sun or the Earth. For the Sun, 
we use the analytic approximations to the results of 0| that can be found in and the solar model we use is that 
of ^3 , with additional abundances taken from . For the Earth, we follow |44j . 

The capture rate of neutralinos inside the Earth receives an additional contribution from a subpopulation of neu- 
tralinos that scatter on a nucleus located near the surface of the Sun, and lose enough energy to stay in Earth-crossing 
orbits which, due to planetary perturbations, do not intersect with the Sun This addition to the local density of 
dark matter in the Earth has a characteristic velocity that more closely matches the escape velocity from the Earth 
than the background halo population, enhancing the resonant capture off elements such as iron. This effect is more 
important for the light, ^ 100 GeV, neutralinos that we are considering and we thus take it into account when 
computing capture rates in the Earth. 

The capture of neutralinos in the Sun or the Earth depends on the elastic scattering cross sections with the nuclei 
that make up the body. These cross sections can be derived j3E| from the nucleon (proton or neutron) cross sections 
that are discussed, for the NMSSM, in App. ^ We have included both spin-independent and spin-dependent terms 
in our computations, the latter being potentially important to evaluate the accretion in the Sun. 

We show in Fig. [21 the equilibrium time between capture and annihilation in the Earth and in the Sun. Most bino 
and mixed-like neutralinos have reached equilibrium and tanh^ {tQ/rA) ~ 1 in Eq. (|ll|l . In the Earth, featuring a 
shallower gravitational potential well, equilibrium has only been reached by a few mixed-like neutralinos, and the 
annihilation rate will be below C/2. The emission region is however much closer to the detector, and, as we will see 
below, and in contrast to the usual situation in the MSSM, the constraints from the Earth are more stringent than 
those from the Sun. 



2. Muon fluxes 

The neutralino annihilation products will hadronize and/or decay giving rise to high energy neutrinos, Ei, < 
m^, which may be detected in a neutrino telescope by measuring the upward-going muons produced by and 
interactions in the rock below the detector. 
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FIG. 2: Equilibrium times between capture and annihilation in the Earth (left) and the Sun (right). Models below the line 
marking the age of the solar system have attained equilibrium and will have an Fann ~ C/2 
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FIG. 3: Integrated muon fluxes above > IGeV from the Earth (left) and the Sun (right). The horizontal line displays the 
MACRO bound \5^. 



The precise determination of the secondary neutrino spectrum is a difficult problem that calls for dedicated Monte 
Carlo simulations of the hadronization and energy losses in the medium of the body. We have used the results of [i^ 
and adapted the relevant routines in DarkSUSY to take into account the additional Higgses present in the 
NMSSM. Spectra are given for six fundamental channels, XX ~^ cc, bb, ti, tt, W^W^ , ZZ, which are also used for 
the Higgs and Higgs/gauge boson final states by following the decay chain until one of the fundamental channels is 
reached. 

The muon yields in include the effects of hadronization/decay of the annihilation products, interactions on 
their way out of the Sun and near the detector, and of the multiple Coulomb scattering of the /i on its way to the 
detector. A similar study was done in [50l |. 

The effects of oscillations in the propagation of the neutrinos through the Sun have been discussed in More 
recently, the full spectra of all neutrino flavors including additional channels, such as light quarks and gluons, and 
accounting for oscillations and ^'t- -regeneration were given in .52]. The combined effect amounts to a 0(0.1 — 10) 
correction which is comparable to astrophysical uncertainties. We do not include these effects here, although if an 
anomalous v signal were discovered, it would be then interesting to try to reconstruct the mass and branching ratios 
of the would-be neutralinos 52] . 

Fig. 13 shows the muon fluxes for xX annihilation from the Earth and from the Sun. We show integrated fluxes 
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above a threshold energy of 1 GeV and the horizontal line represents the MACRO limit |53 which is comparable 
to that of Super-Kamiokande [s^l- To be able to constrain the low mass neutralinos considered in this work, it is 
crucial for the detector to have a low threshold. Of all the forthcoming facilities, ANTARES [s^ seems to be the 
most promising one, with an advertised threshold of i?^ ~ 10 GeV and a target sensitivity of 100 — 1000 km~^yr~-'^, 
it should be able to detect or further constrain those models with > 20 GeV. Larger facilities like ICECUBE [s^ 
have sparser instrumentation, which increases the threshold to i?^ > 100 GeV, above the mass range in which we are 
interested. 

In Fig. 121 we considered a half- aperture 9 — 30°. A cone of this size should contain most of the muons from 
annihilations of even the lightest neutralinos. For moderately larger masses, a smaller aperture could improve the 
signal to noise ratio by reducing the background while still collecting most of the signal, and the limits could be 
improved by an optimized analysis. 

It is encouraging, however, that present muon fluxes due to capture in the Earth, presumably in part due to the 
enhancement in the density of light neutralinos in solar system orbits, are currently able to rule out a few models with 
moderate masses, ~ 60 — 80 GeV, and that an order of magnitude improvement in sensitivity as expected with the 
ANTARES telescope, should enable to access a sizable part of the parameter space by looking at signals from both 
the Sun and the Earth. On the other hand, singlino-like neutralinos show suppressed muon fluxes and prospects for 
their detection seem more remote. 



B. Solar Physics Bounds 

Energy transport by neutralinos could have relevant effects on the Sun, producing an isothermal core and reducing 
the Sun central temperature, Tc- WIMPs with masses of a few GeV and elastic scattering cross-sections around 
aei ~ 10~^^ cm^, were considered some time ago as being able to reduce the solar neutrino flux, hence solving the 
solar neutrino problem Is^ IssL Is^ . It has, since then, been realized that the solar neutrino problem cannot be 
solved by simply reducing Tc, and this hypothesis was abandoned. 

On the other hand, our knowledge of the solar interior has advanced to a point where stellar evolution theory 
in combination with observational data could provide information on the existence and properties of the particles 
constituting the dark matter. The sound speed in the Sun is known with an accuracy of roughly 0.1% through 
helioseismic data 6(1 |, a nd the measurement of the neutrino flux from decay has enabled the determination of Tc 
at the percent level |6ll |. 

The variations in the sound speed induced by dark matter particles were considered in l6l| and, together with 
the influence on the boron neutrino flux 63], were claimed to exclude WIMPs below ^ 60 GeV. This stringent 
conclusions were due, according to (64| . to an unrealistic extrapolation of the helioseismic data down to the central 
regions of the Sun. Neutralinos as light as ~ 30 GeV were shown to be in accord with helioseismology and also to 
leave the neutrino fluxes unchanged, since the central temperature was only being modified in a small region around 
the center of the Sun. 

It is nonetheless of interest to consider the influence on the solar energy transport of neutralinos within the NMSSM. 
Apart from changes in the capture rates, the mass of the neutralinos we are considering here dwell well below 
my, ^ 30 GeV, creating a larger isothermal core with potential observable effects. 

Energy transport in the Sun can occur by diffusion or in a non-local manner. The prevalence of either regime is 
determined by the Knudsen number, which is the ratio of the mean free path of the weakly interacting neutralino in 
the multicomponent baryonic background to the scale length of the system: 

Kn= ^L^2^i'^i^ , (13) 

where the sum runs over the chemical elements in the Sun. 

For neutralinos in the Sun, the relevant geometric dimension is the scale height of the neutralino cloud in the central 
region, which can be approximated by: 



L = r^^0.13i?oji^. (14) 
V "^x 

When the mean free path is short compared to r^, energy is transported by thermal conduction and the relevant 
Boltzmann collision equation has been carefully studied in We will be mostly interested in the opposite regime, 
the Knudsen limit, when Kn 3> 1 and the particles orbit many times in the Sun between interactions with nuclei. 
An analytic approximation for this case was presented in j although Monte Carlo simulations |6^ revealed that it 
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overestimated the neutralino luminosity by a factor of a few. This was subsequently confirmed, and the source of the 
discrepancy attributed to the deviation from isotropy of the neutralino distribution . With this caveat in mind, it 
will suffice, for our purposes, to estimate the neutralino luminosity using the results of [sl ]. 

We asserted above that the neutralinos will be transferring energy in the Knudsen regime: let us show now that 
this is indeed the case. The critical cross section for an interaction to occur in a solar radius can be estimated as: 



36cm2 



(15) 



The x — n,p elastic scattering cross sections that we obtain using the results in Add. IbI fall a few orders of magnitude 
below (Tc- For some mixed-like neutralinos we get values as large as 2 x 10~^^ cm^, an order of magnitude larger than 
for bino-like neutralinos and some three orders of magnitude above those of singlinos. Hence, the neutralinos will 
travel over distances larger than 10"^ Rq, corresponding to Knudsen parameters in the range Kn > 10^. We show the 



parameter in Fig. ^ at a distance from the center of the Sun, which is where the neutralino luminosity is expected 
to peak [67|^°. 




40 60 
m, IGeV] 



100 



FIG. 4: Knudsen number, Eq. at a distance 



The energy transport is most effective in the region Kn ~ 0.4 [6^ |63, far below the values depicted in Fig. ^ 
so we expect the neutralino luminosity to be a fraction of the total solar luminosity Cq. Indeed, we can obtain a 
rough estimate of C^, by adapting Eqs. (2.8-2.10) in js^ to account for the different species of nuclei in the Sun, and 
assuming the neutralino luminosity is confined to a region of size : 



r*Kn 



'iGeV 



4.1 X IO^^Cq, 



(16) 



where the number of neutralinos in the Sun is given by N^^ = CrAtanh(i0/T_A). 

Looking at Fig. |31 where we see that most models contribute a tiny fraction of the solar luminosity. However for 
the lightest bino-like neutralinos, m-^ ^ 30 GeV, the neutralino luminosity may be comparable to the total solar 
luminosity, and may thus already be disfavored. 

We have already mentioned that the neutralino luminosity might be overestimated by a factor of 0(10). On the 
other hand, the neutralino luminosity is not directly observable, and it is not unconceivable that neutralinos giving a 
smallish fraction of the total luminosity, but having a small mass and, hence, a sizable isothermal core, might modify 
the boron neutrino fluxes appreciably. As pointed out in ^5&]. a 20 GeV neutralino carrying only 1O~^>C0, could be 



In this respect, the quantity 5 used in 
luminosity at the center of the Sun. 



IMI does not seem appropriate to characterize neutraUno energy transport, since it involves its 
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FIG. 5: Neutralino luminosity in units of Cq. 



responsible for the transfer of up to 50% of the energy in the inner region bounded by r^. Computing the actual 
modifications in neutrino fluxes and/or helioseismic data, would require the generation of self-consistent solar models 
with the neutralino transport taken into account. Although this is beyond the scope of the present work, the present 
estimates suggest that this task may deserve further investigation for light NMSSM neutralinos. 



C. Antimatter from Neutralino annihilations in the Galactic Halo 



Neutralino pair-annihilations in the galactic halo can produce, through the hadronization or decay of the underlying 
elementary constituents arising from the annihilation process, antimatter in the form of positrons and of hadronic 
stable antimatter states like antiprotons and antideuterons. The abundance of antimatter fluxes produced in neutralino 
pair-annihilations not only depends upon the particle physics nature of neutralinos, but also on various astrophysical 
factors. The latter-including the structure of the dark matter galactic halo, the propagation of cosmic rays in the 
Galaxy, the effects of solar modulation-induce some amount of uncertainty in the flux computation. Further, while in 
the case of low-energy antideuterons the cosmic ray background can be suppressed at a level where the detection of 
even a single antideuteron can be a signal for new physics and potentially for dark matter annihilations in the halo, 
for positrons and antiprotons the background is large. While this latter background is, to some extent, understood, 
it has to be properly incorporated and estimated if one is to be able to extract a possible dark matter annihilation 
signal from the data. 

As far as the dark matter distribution in the galactic halo is concerned, we resort here to the strategy outlined in 
Ref. [H3 (the reader is referred to Ref. jH, [iS ll^ 

for more details). We consider two extreme possibilities for the 
structure of the dark matter halo. In the first scenario, the central cusp in the dark matter halo, as seen in numerical 
simulations, is smoothed out by a significant heating of the cold particles leading to a cored density distribution, 
which has been modeled by the so called Burkert profile |73|] , 

''^'^ = (l + r/a)a + W) - ^''^ 

Here, the length scale parameter has been set to a = 11.7 kpc, while the normalization is adjusted to reproduce 
the local halo density at the Earth position to p_b(^o) — 0.34 GcV cm^"^ 69J. We refer to this model as to the Burkert 
Halo Model. It has been successfully tested against a large sample of rotation curves of spiral galaxies Jt^I . 

In the second scenario we consider here, baryon infall causes a progressive deepening of the gravitational potential 
well at the center of the galaxy, resulting in an increasingly higher concentration of dark matter particles. In the 
circular orbit approximatron © [li , this adiabatic contraction limit has been worked out starting from the N03 
profile proposed in Ref. [77| : the resulting spherical profile, which has no closed analytical form, roughly follows, in 
the inner galactic regions, the behavior of the profile proposed by Moore et al.,^^, approximately scaling as r~^-^ in 
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ANTIDEUTERONS (ULDB Balloon Flight) 
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FIG. 6: The expected number of antideuterons detectable with an ultra-long duration balloon-borne GAPS-type experiment, 
as a function of the lightest neutralino mass. In the left panel we adopt a Burkert halo model, while in the right panel we make 
use of an adiabatically contracted N03 halo profile. The conventions for the various neutralino types follow those of Fig. Q 



the innermost regions, and features a local dark matter density p^vosC^o) — 0.38 GeV cm^'^. We dub this setup as the 
Adiabatically Contracted N03 Halo Model. 

The parameters for both models have been chosen to reproduce a variety of dynamical information, ranging from 
the constraints stemming from the motion of stars in the sun's neighborhood, total mass estimates from the motion of 
the outer satellites, and consistency with the Milky Way rotation curve and measures of the optical depth toward the 
galactic bulge 69, 70]. Both models have been included in the latest public release of the DARKSUSYpackage 

The antimatter yields from neutralino annihilation are then computed following the procedure outlined in Ref. [68| . 
We calculate the neutralino annihilation rates to p and n using the Pythia 6.154 Monte Carlo code [l^ as implemented 
in DarkSUSY and then deduce the D yield using the prescription suggested in Ref. The propagation of 
charged cosmic rays through the galactic magnetic fields is worked out through an effective two-dimensional diffusion 
model in the steady state approximation (Slij , while solar modulation effects were implemented through the analytical 
force- field approximation of Gleeson and Axford 82]. The solar modulation parameter $f is computed from the 
proton cosmic-ray fluxes, and assumed to be charge-independent. The values of <I>i? we make use of refer to a putative 
average of the solar activity over the three years of data taking of the recently launched Payload for Antimatter 
Matter Exploration and Light-nuclei Astrophysics (PAMELA) experiment [s^ for positrons and antiprotons, and 
over the estimated period of data-taking for the General Anti-Particle Spectrometer (GAPS) experiment in the case 
of antideuterons. 

For antideuterons we consider the reach of the proposed general antiparticle spectrometer (GAPS) [s^, El] in an 
ultra-long duration balloon-borne (ULDB) mission, tuned to look for antideuterons in the very low kinetic energy 
interval from 0.1 to 0.25 GeV per nucleon. As described in Ref. [s^, in fact, this experimental setting would allow one to 
safely neglect the background from secondary and tertiary cosmic-ray-produced antideuterons, unlike a satellite-borne 
mission: the detection of a single low-energy antideuteron would then be a clean signature of an exotic antideuteron 
source (including, but not limited to, galactic dark matter annihilation). We set the value of the solar modulation 
parameter at the value corresponding to the projected year for the balloon-borne GAPS mission, around 2011. 
The resulting sensitivity of GAPS has been determined to be of the level of 3 x 10^*m^^sr^^GeV~^s~^ [ssl. Is^. 

To evaluate the sensitivity of the PAMELA antimatter search experiment, we adopt the statistical treatment of 
the antimatter yields introduced in Ref. '68*] (an analogous approach has been proposed for cosmic positron searches 
[S^). Motivated by the fact that the signal is much smaller than the background, we introduce a quantity which 
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FIG. 7: The Visibility Ratio for positrons, as defined in Eq. (I21I I. as a function of the fightest neutralino mass. In the left panel 
we adopt a Burkert halo model, while in the right panel we make use of an adiabatically contracted N03 halo profile. The 
conventions for the various neutralino types are as in Fig. Q 



weighs the signal's "statistical significance, summed over the energy bins" , 



- / „ ^^^^ 



where 4>s{E) and (t>biE) respectively represent the antimatter differential fluxes from neutralino annihilations and from 
the background at antiparticles' kinetic energy E, and Tmin. max correspond to the antiparticle's maximal and minimal 
kinetic energies to which a given experiment is sensitive (in the case of the PAMELA experiment fs^, T^^^^-^ ~ 50 
MeV, T^tx = 270 GeV, T^j^ = 80 MeV and Tf^^^ = 190 GeV). It can be easily verified that Eq. ^ reproduces, in 
the large-number-of-bins limit, the excess from an exotic contribution in the fit to the prospect antimatter fluxes. 
We compute the primary component, 0s, with the DARKSUSYpackage, interfaced with a subroutine implementing 
the diffusion and solar modulation models outlined above. The background flux 0t, has been calculated with the 
Galprop package [S^, with the same propagation and solar modulation parameter choices employed to compute the 
signal. 

Given an experimental facility with a geometrical factor (acceptance) A and a total data-taking time T, it has been 
shown |68j | that, in the limit of a large number of energy bins and of high precision secondary {i.e. background) flux 
determination, a SUSY model giving a primary antimatter flux (j)s can be discriminated at the 95% C.L. if 

/40s).A.T>(x^)^f , (19) 

where (x^)^^" stands for the 95% C.L. with Ub degrees of freedom. For the PAMELA experiment, where A ~ 
24.5 cm^ sr, T=3 years and Ub — 60 we get the following discrimination condition [68j | 

740^) > l^ii^ = /3y- PAMELA, 95% ^ 3 3 ^ lO"* cm-^sr-^s-i (20) 



which is approximately valid for both positrons and antiprotons (though in the latter case the PAMELA experiment 
is expected to do slightly better). As a rule of thumb, the analogous quantity for AMS-02 should improve at least by 
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FIG. 8: The Visibility Ratio for antiprotons, as defined in Eq. 1211 1. as a function of the hghtest neutrahno mass. In the left 
panel we adopt a Burkert halo model, while in the right panel we make use of an adiabatically contracted N03 halo profile. 
The conventions for the various neutrahno types are as in Fig. Q 



one order of magnitude 36]. In our plots, we will show, for both antiprotons and positrons, the following "Visibility 
Ratio" 

Visibility Ratio ^ /f" /if^ ^^^^^^^ ^^^^ . (21) 

We show our results on the prospects for detecting a WIMP pair annihilation signature in the various above 
mentioned antimatter channels in Fig. I6I8I As in the previous figures, we indicate singlino-like models with red 
squares, bino-like with green pluses and mixed singlino-bino models with empty blue circles. Models lying above the 
horizontal lines are expected to give a detectable signature at the experiments discussed above. In all three figures, 
we adopt the Burkert Halo Model in the left panels and the adiabatically contracted N03 Halo Model in the right 
panels. As a general comment, switching from the conservative Burkert profile to the more optimistic adiabatically 
contracted halo profile causes an increase in the fluxes of around one order of magnitude (notice that, in terms of the 
Visibility Ratio, Eq. (|2I|I . for antiprotons and positrons, which depends on the square of the signal flux, this translates 
into a two orders of magnitude increase). 

We start showing, in Fig. |H| the Visibility Ratio for antideuterons, effectively given by the expected number of 
detected antideuterons at an ULDB GAPS mission. As alluded above, this experimental setup is virtually devoid 
of cosmic ray background, hence the detection of even only one D can be regarded as a "signal . We notice that 
in general, low mass neutralinos, peculiar of the NMSSM setup under consideration here, yield a sizable flux of low- 
energy antideuterons. With some exceptions, singlino-like neutralinos produce an insufficient flux of D, while the 
most promising models are mixed singlino-bino models with a mass in the range 55 ^ m^/GeV ^ 95. 

Fig. Hand IHl respectively show the Visibility Ratios for positrons and for antiprotons. As a general comment, we 
point out that in the present setup antiprotons stand as a more promising channel to effectively disentangle an exotic 
signal. As for the case of antideuterons, low mass models are again expected to give a sizable antimatter yield. While 
we do find some instances of singlino-like neutralinos that can give large antimatter fiuxes, in general we find that the 
antiproton and positron yield from singiinos is not particularly promising. On the other hand, mixed models, peculiar 
to the NMSSM, give in general large fiuxes, and a significant portion of the models will be tested by the results from 
the space-based PAMELA experiment on a time scale of three years (or by AMS-02 in a much shorter time scale). 

We also computed the constraints from current antiproton j89| and positron l90| fiux measurements, in terms of 
the to the data of the sum of the background and the signal. Using this criterion, we find that models featuring 
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an antiproton Visibility Ratio larger than ~ 10 are generically conflicting with current data, and so are models 
giving a positron Visibility Ratio larger than ~ 300. However, one should keep in mind that the background we use 
in our computation can be somewhat lowered without conflicting with cosmic ray propagation models; in a more 
conservative approach, asking that the signal alone does not exceed the measured antiproton flux, rules out only 
models with Visibility Ratios larger than ~ 100 (~ 1000 in the case of positrons). 



D. The Monochromatic Gamma- ray Flux 



Neutralinos can pair annihilate in the Galaxy or in dark matter concentrations outside the Galaxy yielding a coherent 
and directional flux of gamma rays; two components add up in the total gamma-ray yield expected from ncutralino 
pair-annihilations: a continuum part, extending up to gamma-ray energies i?^ ^ m^, generated by annihilation 
products radiation and from decays of, e.g., tt*^ —> 77, and (possibly more than one) monochromatic lines, in loop- 
induced direct decays to, e.g., 77, Z7 or final states. Among the latter, the brightest, and the one which occurs in 
any supersymmetric framework (the others being potentially kinematically forbidden) is often that associated to the 
77 final state. Since the possibility of unambiguously disentangling the continuum gamma-ray contribution from the 
background is known to be observationally extremely challenging (see e.g. the recent analyses in Ref. [9lll93 . l9^ 1. and 
in view of our expectations on the size of the 77 annihilation channel in the NMSSM, as anticipated in the Introduction, 
we shall concentrate here on the monochromatic gamma-ray line from radiative annihilation of neutralinos into two 
photons, at an energy E-y — m^. 

As well known, the estimate of the gamma-ray flux from WIMP pair annihilation critically depends upon the 
assumptions one makes on the dark matter profile in the inner portions of the halos. This spread can be extremely large 
in the case of the nearby Galactic Center, where the dark matter distribution is poorly constrained by observational 
data. One is then forced to extrapolate the assumed dark matter profile to very small regions around the center of 
the Galaxy; the small scale central structure of dark matter halos plays, instead, a less crucial role when the source 
is located further away [oJl , as in the case of nearby dwarf satellite galaxies |^ |^, 95] or of galaxy clusters [o^ . 
A second issue involved in the evaluation of the possibility of detecting a WIMP annihilation signal in gamma-ray 
data is related to the evaluation of the background. In short, any evaluation of the detectability of a WIMP induced 
gamma-ray signal must be carefully and properly put in a specific context; comparing the detection perspectives for 
different astrophysical WIMP annihilation locations can be even more difficult, and full details about the assumptions 
involved have to be specified. 

In Fig. Owe compare the detection prospects, in the (m^, {av)^^) plane (where {<Jv)^^ = {av) x BR(xx ~^ 77))i of 
the 77 line from ncutralino pair-annihilations in the Galactic Center (left) and in the Draco dwarf spheroidal galaxy 
(right) with the predictions we obtain in our scan over NMSSM models. We consider the sensitivity of GLAST after 
five years of data takin g ti me T, assuming an average angular sensitivity of Afi ~ 9 x 10~^ sr, and an average effective 
area A^q of 5000 cm^ [93 ■ We consider a putative energy bin centered around the location of the gamma-ray line, 
= m^, and as wide as the expected energy resolution of GLAST, AE/E ~ 0.1. Namely, we consider the energy 
interval 

iAE)m^ = TOx/1.05 <E^< x 1.05. (22) 
Given a background with a differential flux 



dE ^ 1 GeV J 
we obtain, over the considered energy range a total background flux of 



^ 7 - 1 \ V 1 GeV J V 1 GeV 

The signal flux from the monochromatic line is instead given by 



(23) 



cm-^s-^ (24) 



- ^ ( lO-^W.-. 'l ( -Hl'.m - An ™-V'. (25) 



In the formula above, we defined the dimensionless quantity 

J(^,A1])^-1- / d»— j— f ^ 3 ) / Plumm- (26) 

Anj^^ 8.5 kpc V0.3 GeV/cm-iy Jiinc of sight 
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FIG. 9; Prospects for the detection, with GLAST, of the monochromatic 77 gamma-ray line in the Galactic Center region 
(left) and in the Draco dSph (right). See the text for details. The conventions for the various neutralino types are as in Fig.0 



We define a signal as "detectoWe" provided the number of signal events in the considered energy bin is larger than 
5, and the foUov^ring 5 — tr significance condition is fulfilled: 

Evaluating the gamma-ray background in the Galactic Center is certainly a non trivial task. Since the EGRET data 
from the Galactic Center likely include a gamma-ray source with a significant offset with respect to the actual Galactic 
Center [o^, we shall consider here the data from the HESS collaboration [o^, which feature a much better angular 
resolution. The HESS data from the Galactic Center region indicate a steady power law gamma-ray source with a 
spectrum dN^/dE oc E~^ '^ extending over a range of gamma-ray energies of almost two orders of magnitude |99j| . 
The flux at low energies, ~ 200 GeV, is limited by the experimental energy threshold. Extrapolating down to the 
energies of interest here (a few GeV up to 100 GeV) involves invoking a particular nature for the mechanism responsible 
for the gamma-ray production. Following "931, we consider two extreme choices for the background extrapolation at 
lower energies, namely the models number 2 and 3 of Aharonian and Neronov, Ref. [lOQj . (we shall indicate hereafter 
the two models as A-N2 and A-N3), respectively giving the smallest and the largest extrapolated background levels 
among those considered in |9^. Model A-N2 invokes inelastic proton-proton collisions of multi-TeV protons in 
the central super-massive black-hole accretion disk, while model A-N3 results from curvature and inverse Compton 
radiation. We assume 0o — 1 x 10" cm s" ^^GeV"^ and 7 = 2.0 for model A-N2, while 0o - 3 x 10""^ cm^s ^GeV 
and 7 = 2.75 for model A-N3. 

As far as the values of J(0, Afi) are concerned, we consider the range given by the extrapolation of the two halo 
models considered above (the Burkert and the adiabatically contracted N03 profiles), giving, roughly, J ~ 10 and 
J ~ 10^. The left panel of Fig. |51 illustrates our results. All the sensitivity lines correspond to the criterion given in 
(|27|l . which we find to be always more stringent than the > 5 requirement. The solid blue lines correspond to the 
A-N2 background model, while the A-N3 background is assumed for the red dashed lines. Our results show that even 
assuming a very optimistic dark matter profile, the likelihood of obtaining a significant gamma-ray line detection from 
the Galactic Center is rather low; at the best, a weak excess can be detected either with very low mass neutralinos, 
or with mixed neutralino-singlino models with a mass ~ 60 80 GeV. 

In the case of Draco, the estimated background is only given by the diffuse gamma-ray back grou nd, which we 
parameterize with (j)o ~ 6.3 x lO"^'^ cm^s~^GeV~^ and 7 — 2.1. We follow the results of Ref. [23] as far as the 
estimate of J are concerned; conservatively, a range of viable halo profiles for Draco gives 10 < J < 100. Taking into 
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FIG. 10: Prospects for the detection, with GLAST, of the monochromatic 77 gamma-ray line in the extra-galactic gamma-ray 
background. See the text for details. The conventions for the various neutralino types are as in Fig. 



account the possibility of a central super-massive black-hole and the subsequent adiabatic of dark matter in a central 
"spzfce" can greatly enhance the viable values of J, up to the level of 10^, a value we assume for the black solid 
line. As for the Galactic Center, the prospects of cleanly detecting a gamma-ray line from the direction of Draco do 
not seem particularly exciting, although, again, some models might in principle, and very optimistically, give some 
evidence of an energy-localized gamma-ray excess. 

The monochromatic WIMP pair annihilation is also constrained by the contribution that ann ihilations occurring 
in any dark matter halo and at all redshifts g ive t o the extragalactic gamma-ray radiation |lOll We refer the 

reader to the thorough discussion given in Ref. |l02j . and we make use here of the constraints, on the m^, {(tv)^^ plane 
derived in Fig. 15 of the same study. In particular, we report in Fig. 1101 left, the sensitivity, on the above mentioned 
plane, expected from GL AST , under the two extreme scenarios for the halo profiles and t he p resence of dark matter 
substructures outlined in |102| . The upper curve refers to halos modeled by a NFW profile |l03j | , no substructures and 
concentration parameters inferred from the Bullock et al. model [iQ^, while the lower curve assumes the (cuspier) 
Moore et al. profile [tH, with 5% of the halo mass in substructures with concentration parametes 4 times than that 
estimated with the BuUok et al. model. In the most generous scenario, a few mixed singlino-bino models can give rise 
to a detectable signal at GLAST, although more conservative assumptions leave small space for any hope of detecting 
any signature at all in the extra-galactic gamma-ray data. 

Even though the prospects for the detection of the monochromatic line do not look particularly promising here, we 
wish to point out that the branching fractions we find, and the absolute values of {av)^^ are, typically, larger than in 
the MSSM. We devote Add. I A 2"al to a detailed discussion of this point, but we wish here to emphasize the main reason 
why the NMSSM rate for the process XX ~^ 11 is expected to be more significant than in the MSSM. The potentially 
light extra CP-odd gauge boson gives rise to the extra contributions shown in Fig. ^2 the size of this contribution, 
generically, depends upon whether the annihilation proceeds close to the s-channel resonance (m^ ~ TOaj/2). We 
illustrate the effect of the extra NMSSM diagrams in Fig. ^| right, where we show the size of BR(xx 77) as a 
function of the ratio nia^l {^m^). As evident from the figure, the largest branching ratios occur when rriaj^ / {2m^) ~ 1, 
and they are more than a couple of orders of magnitude larger than the MSSM limit (bino-like neutralinos and 
maj{2m^) > 1). 
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IV. CONCLUSION 

Gauge singlet extensions of the Higgs sector of the minimal supersymmetric Standard Model provide well motivated 
theoretical and phcnomcnological laboratories. Besides offering an elegant solution to the supersymmetric /i problem, 
they provide a viable way out of the difficulties connected to encompassing a mechanism of electroweak baryogenesis in 
the MSSM. In the present analysis, we focused on one specific such extension, the so-called NMSSM, and investigated, 
for the first time, the prospects for neutralino dark matter indirect detection. 

The present study is motivated by two basic observations: first, in the NMSSM, unlike the MSSM, the lightest 
neutralino can be naturally very light, as a result of the possibility of it annihilating through a potentially very light 
extra CP-odd, mostly singlet-like Higgs boson; second, the extended Higgs sector leads to extra diagrams in the 
loop-amplitude relevant for the pair annihilation of neutralinos in photon or gluon pairs. An enhancement of the 
monochromatic XX ~^ 77 gamma-ray line is therefore generically expected within the NMSSM, as opposed to the 
minimal supersymmetric setup. 

We found that the rate of neutrinos produced by the annihilation of neutralino dark matter particles captured 
inside the Earth and the Sun is in general large in the NMSSM; unlike the MSSM, we found that most models give 
a larger signal from annihilations in the core of the Earth rather than in the Sun, at a level which can, in certain 
cases, be constrained by current available data from SuperKamiokande and MACRO. This is presumably due in part 
to additional low velocity contributions to the local neutralino density in the region of the Earth which can result for 
light neutralinos. Future neutrino telescopes with increased sensitivity for low neutrino energies will be able to probe 
a sizable part of the parameter space by looking at signals from both the Earth and the Sun. 

The dynamics of the Sun could also be modified due to energy transport by neutralinos. Our estimates show 
that, especially bino-like, neutralinos below < 30 GeV might contribute a significant fraction of the total solar 
luminosity. More detailed studies, using self-consistent solar models, could unveil large enough modifications on the 
sound speed or on the boron neutrino flux to significantly disfavor light neutralino scenarios. 

We showed that within the NMSSM the expected antimatter yield from neutralino pair annihilations in the galactic 
halo can be sizable, although the absolute normalization of the flux depends on speciflc assumptions about the dark 
matter halo profile. In particular, we found that signals at low-energy antideuteron search experiments such as GAPS, 
and at space-based antimatter search experiments such as PAMELA, are expected, though not guaranteed, for very 
light neutralinos (m^ ^ 20 GcV) or for intermediate mass mixed singlino-bino neutralinos (60 ^ ^ 90 GeV). 

Wc worked out for the first time the loop-induced pair annihilation cross section for NMSSM neutralinos into two 
photons and two gluons, pointing out that the expected branching ratio, with respect to tree-level neutralino pair 
annihilation into other Standard Model particles, is typically large, especially when compared to the MSSM case. 
The reason for this enhancement is traced back to diagrams which are resonant when ~ , the latter quantity 
indicating the mass of the lightest, extra CP-odd Higgs boson. 

We analyzed in detail the prospects for the detection of the monochromatic gamma-ray line resulting from xX ~^ 11 
annihilation processes in the Galactic Center, in a nearby dwarf spheroidal Galaxy (Draco) and the coherent effect 
of annihilations in any dark matter halo contributing to the extra-galactic gamma-ray radiation. Wc pointed out 
that most models are not expected to give any detectable signal at GLAST, although this detection channel looks 
significantly more promising than in the usual MSSM setup. 

Finally, with the purpose of making the present study a useful and complete starting point for future research in the 
field, and in order to sort out and clarify some notational ambiguities and inconsistencies, we collect in the Appendix 
the details of the one-loop computation of the xX ~^ 11: 99 amplitudes and other quantities relevant for the estimate 
of indirect detection rates. 



Acknowledgments 

We gratefully acknowledge useful conversations with John Beacom, Marco Cirelli, Ulrich Ellwanger, Cyril Hugonie, 
Bob McElrath, Alexander Pukhov and Miguel Angel Sanchis-Lozano. 

FF and LMK are supported in part by grants from the DOE and NSF at Case Western Reserve University. SP is 
supported in part by DOE grants DE-FG03-92-ER40701 and FG02-05ER41361 and NASA NNG05GF69G. 



19 




FIG. 11: Additional diagrams, due to the second CP-odd Higgs boson, for the process XiXi ~^ 77 in the NMSSM. The 
contribution of ai corresponds to Fig. l.d (left) and 2.d (right) in ,2^ for the MSSM. 



APPENDIX A: NEUTRALINO ANNIHILATION CHANNELS 
1. Tree level processes 

Analytic expressions for the NMSSM-like neutralino annihilation into two particles at tree level can be found 
in [3". 

For the study of indirect detection, only the nonzero terms in the limit u — s- need to be taken into account, which 
restricts the relevant processes to those with a CP-odd final state that have a non-vanishing S'-wave amplitude. 
With respect to the MSSM, the main differences are: 

• An additional scalar Higgs, h^, exchange in the s channel contributes to the annihilation to W^W— and ZZ. 

• A fifth neutralino is exchanged in the t and u channels for the ZZ process. 

• There is an extra Z - scalar Higgs final state: Zh^. Two Higgs pseudo-scalars instead of one, and five neutralinos 
contribute to these reactions. 

• For the W^H^ final state, one has to take into account the contribution of the extra and a2 Higgses. 

• There are five additional final states with a scalar and a pseudoscalar Higgs. Diagrams Xg have to be considered. 

• Finally, for the / — / final state, one has to include the exchange of the additional Higgses and 02. 

On top of that, the different couplings have contributions proportional to A and n not in the MSSM 13]. We used 
the w-independent part of the S'-wave terms in |l6j | for our predictions in Section lllll 



2. One loop processes 

a. Neutralino annihilation into two photons 

Some radiative processes, even if loop suppressed, are of interest for dark matter detection. The annihilation to two 
photons, XiXi ~^ 111 has a characteristic monochromatic signature at = m^o/2. This allows a clear distinction 
from all astrophysical backgrounds, unlike the continuum spectrum produced i n tree level processes. 

In the context of the MSSM, a full one-loop calculation was performed in p5ll05l |. We computed the cross-section 
for this process in the NMSSM by adapting the results of '23|^^. 

Apart from the dependence on A and n of the NMSSM couplings, we need to compute two additional diagrams, 
shown in Fig. ^2 due to the presence of a second pseudoscalar Higgs boson, 02. 

Four types of Feynman diagrams, contributing to the two photon annihilation amplitude, were identified in [26j |. 
Let us discuss their computation in the NMSSM in turn: 



The notation in lit! follows the usual MSSM practice of labeHng the CP-even Higgs mass eigenstates as H\ and H2- To make contact 
with our conventions one has to switch the indices 1 2 in the scalar Higgs matrix Sij, and 3 <-» 4 in the neutralino matrix Nij. 
Furthermore, A and k in the superpotential have the opposite sign as ours^ 

The expressions for the MSSM were implemented in the code DarkSUSY |49| . We have adapted the relevant subroutines to Nmhdecay. 
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a. Diagrams l.a-l.d For the fermion-sfermion loop diagrams, we need to duplicate the CP-odd Higgs terms 
(Fig.inj and substitute the NMSSM couphngs in Sfj, Dfj, Gz/ and G^,/ of Eqs. (7) and (8) in p6|. 
For up-type quarks, let us define: 



gii 



9rl 



-ff2"lqiVi3 

\/2'mw sin fi 
g2Ni2 + (±2e, - l)gyNii/3 



V2 

gir = V^gyNiiCq 



[Jrr 



an (Al) 



where 52 is the electroweak coupling constant, gy — g2teLii0w, &q is the quark charge and mq its mass. For down 
quarks and Icptons we need to replace Ni^ — > and sin/3 —> cos/3 and use the lower sign for g^i- Then, with: 

91 = cos 6*5 4- .gir sin 6*5 

.92 = -grl sin 9q + grr COS 9q (A2) 

where q is the squark mixing angle which is taken to be q = {i.e. no mixing) for the first two families in Nmhdecay. 
we have: 

f 



2 



Dff = gig2 (A3) 

The ZxiXi coupling reads: 

where the weak isospin, T^f is -1-1/2 for up -typ e quarks and —1/2 for down quarks and leptons. 
Finally, using the aiXiXi coupling from |13l |: 

• ^ V2 

+ ^f{PalIllf - P^^nlf) - ^{PalU^^ - Pa2n^f), (A5) 

where Hfj' = NiaNjb + N.bNja, we obtain: 

^ „ mqg2Pii , , 

^i^iTOu,Sm/5 

Changing sin/3 — > cos/3 and P^i — > Pj2 in Eq. ljA6|l . leads to the corresponding expression for down-type quarks and 
leptons. 

h. Diagrams 2.a-2.d For the chargino-Higgs loop diagrams need also to take into account the additional contri- 
bution of 02 and use the expressions below for S^h, D^h, Gzx ^^"^ ^oix (9) Ell- 
Taking the H~^XiX~ couphng from jl^ : 

sin 3 

dH+xTx" = ^ PUi2Njr-, —U,2igyNji + g2Nj2) + .92 sin/3C/aA^j4 



V2 
cos/3, 

71" 



^ff-x+x" = >^smf3V^2Nj5 + ^fv^2igyNji + g2Nj2) + g2COsf3VaNj3, (A7) 



where [/ and V are the chargino mass matrices, we get: 

^H+x-xl ' ^H-xtx°, 



S 



2 

^xff = 3i/+x-x?5if-x+x?- (A8) 
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The Z exchange diagrams require: 



.y2 / ta2 I tt2 '^i2 \ / Ar2 Ar2 



Gz^ - ( ^/i + ^ - C/A - ^ ) (A^i^3 - Nt,) (A9) 



With Eq. IjASp and the UaxtXj couphng: 



Vx+X- = ^^a3t/^2F,2 - -^(P,iC/,iy,2 + Pa2U,2Vjl), (AlO) 



we have: 



c. Diagrams 3.a-3.c The M/"'"x5x7 couphngs are: 

9wu = 92[~Nl3Va/V2 + Nl2V^2) 

(NuUa/^+Ni2U,2) (A12) 



R 

9wii = 52 



which we can substitute in 5*;,,^^ ^nd D^iy of Eqs. (11) and (12) in ||2^] to obtain the chargino-W loop contribution. 

d. Diagrams 4-o.-4-b The unphysical Higgs boson is orthogonal to the charged Higgs and we can derive its required 
coupling to neutralinos and charginos by adapting those of the charged Higgs in [iJI : 

cos [3 

9g+x-x° " Asin/3f/j2^j5 -^Ui2{9yNji + g2N.i2) + 92 cos f3UiiNj4 

sill 3 

9g~x+x" = Xcos(3VaNj5^—^V^2{9yNji^g2Nj2)+g2Siril5VaNj3. (A13) 



V2 



Then, in Eq. (13) of [2g we need to input: 



SxG - 



DxG = dG^x-xldG-xtxr (^14) 

In the expressions above, we have not taken into account that Nmhdecay uses a real neutralino and chargino 
mass matrix, whereas the expressions in 2(>| assume that diagonalization in the neutralino and chargino sectors is 
performed using a complex iV, U and V , so that m^o and the chargino masses are always positive. 

To correct for this fact we need to multiply by e = sign {rn^o^ all instances of N* in j26j for a vertex in which 

the neutralino is annihila ted [ 2^ and a similar change of sign needs to be done in V . Details for each vertex can be 
found, for the MSSM, in |l06l |. For the two photon amplitude computation, this prescription amounts to multiplying 
by e (77 = sign(m^om^+)) the D terms in diagrams of type 1 (2, 3 and 4). 

The presence of extra light CP-odd Higgses, can enhance the cross section for this process. In Fig. H12|l . we show 
the branching ratio for the process x\x\ ~* 77 together with the contribution from the diagrams in Fig. 1111 

The branching ratio peaks for neutralino masses 2m^o ~ mai, more so, for singiino and mixed-like neutralinos^^. 
The contribution of the CP-odd Higgs diagrams. Fig. II II to the total Feynman amplitude, A, is displayed in the lower 
panels by the quantity 

A ^ ^- — . (A15) 



Larger values of A, corresponding to larger relative contributions of the CP-odd diagrams to the total amplitude, also 
cluster for those values where the branching ratio is larger. 

The light CP-odd Higgses, together with the additional singiino component, lead, thus, to an enhancement of the 
77 annihilation channel in the NMSSM compared to the MSSM. 



Note the absence of singlino-like neutralinos with rn~a < ma-i/'i. 
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FIG. 12: Branching ratio for the annihilation to 77 (up) and contribution of the CP-odd Higgs exchange diagrams to the real 
(left) and imaginary (right) part of the amplitude. 



6. Annihilation into two gluons 

The cross section for this process |l07j can be obtained at once from the two photon channel computed in the 
previous section. In order to do so, we need to consider only the diagrams of type 1 for quarks, with no contribution 
from leptons. The electric charge is substituted by — s- 1, and the color sum average is performed by — > 2a^ in 
the final expression for av |26|. 



APPENDIX B: ELASTIC SCATTERING CROSS SECTIONS 



We review here the computation of the neutralino-nucleon elastic cross-section, which is used to predict direct 
detection rates and, in the context of indirect detection, determines the capture rate of neutralinos in the Sun or in 
the Earth. 

The basic ingredient for the neutralino-nucleon cross-section is the individual neutralino-quark cross-section, which 
for the MSSM can be found in |108J. The x -~ 9 process has been studied before in the context of the NMSSM. First 
in [20|. where both spin- independent and spin-dependent contributions were computed. More recently the problem 
was revisited in [23|, were only the spin- independent part was considered and a mistake in |23| was corrected. The 
authors in j2^ approximated the spin-independent interaction by assuming that the ^-channel exchange of CP-even 
Higgses dominates. For our predictions in Sec. lIIII we re-derived the Xi — g cross-section for the NMSSM by extending 
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the MSSM calculation '105 

The low energy x — g effective lagrangian can be written as: 

^eff = dqxi''ibxqiiiibq + fqxxm (Bi) 

where only contributions that don't vanish when u — > have been written. The first term describes the spin-dependent 
contribution and the second one the spin-independent one. 

As in the MSSM we have two types of diagrams contributing to the spin-dependent interaction in the NMSSM, 
Z exchange and squark exchange. For the spin-independent one, CP-even Higgs exchange and squark exchange 
contribute t o /„. 

Following |108| |. let us define: 

X = -V2[g2T3fNi2-gy{T3f-eq)Nn] 
Y = V2gyeqNn 

_ _Jj2mqNi3_ 7 _ ffamgA^u 
V2sinpm„ cos pm„ 



Then the couplings involving the lightest squark can be written as: 

1 



2 [cos eg{X + Z) + sin 0q{Y + Z)] 



= -[cosdgiX ~ Z) + sin9g{Z -Y)], (B3) 



and the corresponding equations for the heavier eigenstate, j = 2, are found taking sin9q cos 9q and cos 6q — > 
— sin 9q . 

With that, the spin-dependent Xi^l interaction is given by: 

- i E . I^^i: ,2 - ^T^ql (A^?3 " Nl) , (B4) 
- i'^x + '^q) ^ 

where the sum runs over the squark eigenstates and the last term describes the Z exchange contribution. 
As for the spin-independent interaction, we have: 



■^9 ^"7.^ 2 ^ 7 -2-™9 2^ 2 rn^a\-nR- ^^^^ 

* ]~i ~ i'^x + '^q) J~i ^Hj sm p 

Note that in the NMSSM we have three CP-even Higgses, included in the last term. For down type quarks, we need 
to replace in Eq. I)B5|1 . Sji — > 5*^2 and sin/3 — > cos/?. Also, we need the coupling q — q — Xi- 

9h^X°X° ~ + Sa2^ij + SaS^if) ^ y^l-^SasNi^Nj^ 



_ -V "V I — \ 



-^{^al^if - Sa2'nlf) + —{SalUff - Sa2^ij)- (B6) 



Our expression for the spin-dependent interaction agrees with the computation in pol | . As for the spin-independent 
part, the authors in noted a mistake in the expressions given in |23|. We agree with their remark, but note that 
the q — q — Xi couplings given in |2lj | contain a mistake, since the couplings for q2 cannot be obtained from those of 
Qi by the usual change sin 9q cos 9q and cos 9q — sin 9q. Indeed, the sign in front of N^2 should affect the whole 
coefficient of the sin^, term in their Eq. (A. 11). 

Since we have used a real neutralino matrix Nij, we should add the necessary factors of e in the expressions above. 
We have followed 108], where a real iV^ was used, and where it was pointed out that the absolute value should be 
used in the kinematic factors appearing in various denominators. However, following the prescription in p^ . one can 
check that when m^o < 0, one should also take fq efq. This extra sign difference, however, does not affect the 
nucleon cross sections discussed below, for they depend quadratically on fq or dq. 

Once the individual Xi~q cross sections are determined, we can proceed to compute the nucleon (proton or neutron) 
cross sections used in Sec. IIIII 
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The spin-independent nucleon-neutralino elastic cross section is given by: 

4m2 



a. 



V r2 



TT 

where 



/;,„ (B7) 



E + ^ 1- E E frp.^ (B8) 

q—u,d,s ^ \ q—u,d,s j q—c^b,t ^ 

and for the quark composition of each nucleon, f^'^ , we use the central values found in [s^- In Eq. (|B7|I . the reduced 
mass is = mp,nm~^/ {mp,n + "^x)- 

For the axial-vector interactions we need the nucleon spin carried by each quark. We use, again, the central values 
from (3Jj to find: 

< = —3 (/«AS'" + h^T + fs^rf ■ (B9) 

TT 

In Sec, nil B1 scalar cross sections with nuclei are used. Since the values at zero momentum transfer are good enough 
for the task at hand, we compute them as: 



4m: 



2 



(J,. 



'^-[ZU + {A-Z)U]\ (BIO) 



TT 



where Z and A are the atomic and mass numbers of the nucleus. 
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